Tomographic femtosecond X-ray diffractive imaging 
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A method is proposed for obtaining three simultaneous projections of a target from a single 
radiation pulse, which also allows the relative orientation of successive targets to be determined. 
The method has application to femtosecond X-ray diffraction, and does not require solution of the 
phase problem. We show that the principal axes of a compact charge-density distribution can be 
obtained from projections of its autocorrelation function, which is directly accessible in diffraction 
experiments. The results may have more general application to time resolved tomographic pump- 
probe experiments and time-series imaging. 



Unless crystallographic redundancy can be taken ad- 
vantage of, radiation damage provides a well established 
limit to resolution for imaging in biology. X-ray mi- 
croscopy of proteins is limited in this way to 10 nm[l[. 
Since dose depends on the inverse fourth power of reso- 
lution, a severe penalty attends any attempt to improve 
resolution beyond this barrier, which occurs when the re- 
quired dose needed to distinguish adjacent image voxels 
with statistical significance exceeds the damage limit at 
that resolution (voxel size). 

It has been suggested that the development of the free- 
electron X-ray laser (FEL) may break this nexus between 
dose and resolution [2[ , if it can provide sufficient photons 
for a useful diffraction pattern in a single pulse, which 
terminates prior to any of the characteristic times for 
damage processes. A continuum of such times, femtosec- 
ond for electrons, hundreds of femtoseconds for nuclear 
motion, is associated with the various irreversible damage 
mechanisms and excitations. Theory suggests that pulses 
shorter than the Auger decay time of a few femtoseconds 
may terminate before significant disruption of the valence 
electron distribution occurs. Experimental evidence, us- 
ing a 25 fs pulse of soft (30 nm) X-rays, now exists for 
this process of "diffraction-before-destruction" at low (90 
nm) Q and higher periodically averaged 0] resolution. 
Since the FEL generates in excess of 10 12 fully coherent 
photons in such a pulse, the method of diffractive lensless 
imaging [5] , in which real-space images are reconstructed 
computationally from these scattering patterns, would 
appear to provide a means of overcoming the radiation 
damage barrier to high resolution imaging in biology [6| . 
The shortest FEL wavelength is currently 7 nm at the 
Hamburg FLASH facility. Shorter wavelength FELs are 
being planned and constructed at other sites around the 
world. 

The subsequent destruction of the sample following the 
initial elastic scattering event, however, has precluded the 
possibility of three-dimensional (tomographic) imaging of 
unique structures. To overcome this limitation, we sug- 
gest here a means to determine the relative orientation 



FIG. 1: Scheme for tomographic femtosecond diffraction, 
drawn for only two beams for simplicity. Beamsplitter XI is 
set to the dynamical 3-beam diffraction condition. Crystals 
X2 and X3 operate at the 2-beam dynamical condition. KBI 
and KB2 are focusing optics for the target at B, with area 
detectors CCD I and CCD2. 



of successive targets. 

Two possible arrangements are shown in Figs. 1 and 2. 
In Fig. [TJ a beam-splitter and reflecting crystals direct 
three orthogonal beams onto a non-periodic target par- 
ticle producing three far-field diffraction patterns prior 
to destruction of the target. All three two-dimensional 
patterns are read out after each X-ray pulse, whereupon 
a new, identical target such as a biomolecule is inserted 
in a new orientation. 

We will show that the relative orientation of successive 
targets can be determined even if the structure of the 
target is unknown. 

Following related work by Smilgies on crystals 0, we 
propose experimental determination of the principal axes 
of the molecule to describe its orientation relative to the 
laboratory frame defined by the incident probe beams. 
Multiple scattering will be neglected, i.e., the first Born 
approximation is assumed valid, so that the patterns have 
inversion symmetry and the target density is a real func- 
tion. 

To expose the principle of the method, assume that 
the phase problem can be solved, by, for example, it- 
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FIG. 2: Two-beam beamsplitter with sample shown at P 
lying on the exit face of the beamsplitter. The source S is fo- 
cused onto two area detectors Dl and D2 containing central 
beam-dump holes. The two vertical arrows show the direc- 
tion of the Poynting vector. Three such orthogonal diffracted 
beams, rather than the two shown, are proposed in the text. 



erative methods (see Q for a review). (We will relax 
this assumption later.) Then, in the "flat" Ewald sphere 
approximation, i.e. at high probe energy, each beam de- 
livers a projection in real space (along the direction of the 
corresponding beam) of the scattering strength per unit 
volume. The projections will be referred to different (ran- 
domly positioned) origins, and both enantiomorphs (re- 
lated by inversion symmetry) will be present with equal 
likelihood. However once a particular enantiomorph is 
chosen for one projection, the resulting two-dimensional 
envelope will constrain the choice of enantiomorph for 
the other two projections. 

Consider the moments of the mass density p(r) for the 
target Q. The zeroth moment delivers the total mass, the 
first moment delivers the center of mass vector, and the 
second moment delivers the moment of inertia tensor. By 
diagonalising the latter, the principal axes of the target 
may be found and hence its orientation relative to the lab 
frame. Taking the center of mass position as the origin, 
the inertia tensor is 



P{r){ 
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where E is the unit tensor and rr is the outer product 
of the position vector with itself. As with any symmetric 
tensor, I has only six independent elements, real eigen- 
values, and orthogonal eigenvectors corresponding to dif- 
ferent eigenvalues. 

Now take p(r) to be the electronic charge density of the 
target, whose projections in three orthogonal directions 
are provided by the phased data, and which define the 
x, y, and z directions specified by unit vectors e, in the 
lab reference frame. The six independent elements of this 



charge density "inertia" tensor then have the form 

I xy = - Jj p z (x,y)xy dxdy (2) 

and similarly for I xx , I yy , I yz and I xz . Here p a is the 
projected density along the cv-direction. Two of these 
six tensor elements can be computed from each of the 
three projections, e.g. I zz and I xy from the projection 
along the z-axis. Hence the inertia tensor of the target 
is fully specified by computing moments and products of 
inertia from the three projections. This charge density 
"inertia" tensor will differ from that based on mass but 
serves equally well to provide a consistent set of principal 
axes fixed to the molecule and defining alignment. Being 
symmetric, this tensor may be diagonalized by solving 
the eigenvalue equations 



I B — bB 



(3) 



for the three eigenvalues b and corresponding eigenvectors 
B. These eigenvectors define a new orthogonal coordi- 
nate system e^- in which the three unit vectors lie along 
the principal axes of the inertia tensor. Barring degen- 
eracy among the eigenvalues, the three eigenvectors are 
unique to within a sign, and therefore offer a natural 
means of specifying the orientation of the target relative 
to the incident beam directions (lab frame). With the 
unit vectors ej of the lab frame and e^- both known, the 
angles between the principal axes of the target and the 
lab frame can immediately be computed. Thus the ori- 
entation of the target has not only been defined by intro- 
ducing the principal axes of the inertia tensor, but also 
specified (within polarity) relative to the lab coordinates. 

Clearly, then, to establish the orientation of the par- 
ticle it suffices to (i) Record three diffraction patterns, 
one for each of the three incident beam directions, (ii) 
fnvert the diffraction patterns using phase retrieval tech- 
niques to yield three real-space projections of the scat- 
tering strength, (iii) Compute the first moment of each 
projection to obtain the center of mass position for that 
projection, (iv) Compute the second order moments of 
each projection (products of inertia) about the center of 
mass to obtain one diagonal and one off-diagonal tensor 
element, (v) Diagonalize the resulting tensor to obtain 
the eigenvectors of the tensor, (vi) Compute the orien- 
tation of each beam relative to the eigenvectors of the 
target in order to determine the angles between labora- 
tory and principal axes coordinates, (vii) If this process is 
repeated for many successive identical targets in random 
orientations, their relative orientations can be found, and 
hence a complete three-dimensional tomographic image 
can be assembled by standard tomographic techniques 
such as filtered backprojection. 

We now extend this analysis to show that the principal 
axes may be found, even without solving the phase prob- 
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lem, by working with the autocorrelation of the sample 
density 



A{r) 



dr' p(r + r')p(r') 



(4) 



A typical product of inertia is 
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= J dr'p(r') [l p xy - x'y'M] = 2M/^(5) 




where M = J drp(r), and with use of the parallel axis 
theorem to calculate the product of inertia for the shifted 
coordinates. Therefore, the principal axes of the autocor- 
relation function are the same as the principal axes of the 
corresponding density. 

Given a "flat" Ewald sphere, the Fourier transform 
of each diffraction pattern (intensity) directly provides a 
projection of the three-dimensional autocorrelation func- 
tion of the density, and the analysis simply requires 
changing p(r) to A(r) in Eqs. [H 

More generally, the moment of inertia can be calcu- 
lated from the second derivative of the Fourier transform. 
Denoting these by a tilde, 

p(q) = J dre-^ r p(r) , A{q) = p{q)p{-q) 



l£u = dq.dqyMl) _ = d qx d qy p(q)p(-q)\ 
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as before. Replacing A(q x , q y , 0), ("flat" Ewald sphere) 
with the correct Ewald sphere (finite radius) diffraction 
pattern, and taking the incident wave vector k along z, 



dq.dqyA (qx,q y , \Jk 2 -q 2 x -ql - k 
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so that the moments of the Fourier transform of the 
diffracted data still yield the same principal axes. 

We have investigated this procedure using detailed nu- 
merical simulations based on data in the Protein Data 
base for GroEL-GroES protein complex (PDB entry 
1SVT). The three-dimensional density was synthesised 
from the tabulated atomic coordinates. Fig. [3] shows 
the projected densities and corresponding projected au- 
tocorrelation functions using the principal axes obtained 
from Eqs. O A second density was then generated in 
a random orientation with respect to the first, as shown 
in Fig. [4j For each of these orientations the principal 
axes were determined using both the densities and the 
autocorrelation functions, giving similar results. The ro- 
tation matrix needed to rotate from the first (Fig. [3]) 
to the second (Fig. orientation was generated from 



FIG. 3: The three orthogonal projections of the GroEL 
charge density (upper) and the corresponding projections of 
the autocorrelation function (lower). The bar indicates 10 
nm. 






FIG. 4: Projections of GroEL density (upper) and autocor- 
relation function (lower) in a second random orientation. 



the principal axes. As a result of inversion symmetry 
there are four distinct choices of rotation matrix (corre- 
sponding to choices of eigenvector signs) when the auto- 
correlation function is used. The correct rotation matrix 
was obtained by testing each to see which predicted lines 
of intensity in diffraction patterns common to two ori- 
entations. (Any two planes in reciprocal space passing 
through the origin must intersect along a common line). 
In this way only one rotation matrix will be found to 
give consistent results. Numerical trials have found this 
procedure to be reliable with several different test ob- 
jects. Our use of specific common lines should be more 
robust for noisy data than common line search methods. 
(Shncerson et al. [l0( have shown that identifications 
of common lines in diffraction pattern down to a mean 
photon count of 10 per pixel enables the determination of 
their relative orientations without the need to solve the 
phase problem.) 

This treatment easily extends to the case where the 
beams from the beamsplitter are not orthogonal. Re- 
ciprocal vectors can be defined in the usual way, so that 
each pattern lies in the plane of two of these vectors. The 
products of inertia may be simply evaluated in terms of 
these reciprocal vectors, and finally transformed into the 
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required lab frame moments. 

Two separate experimental implementations of this ap- 
proach are suggested in Figs. [T] and [21 each with an in- 
cident beam close to the [111] direction of a diamond 
beam-splitter crystal, set to simultaneously excite the 
[022] and [220] reflections. For an X-ray energy close to 
7 keV, this generates three orthogonal beams with Bragg 
angles of 45°. The Borrmann effect [nj may then be 
used to produce three beams of approximately equal in- 
tensity for a crystal thickness of order one mm, due to 
3-beam multiple scattering 12, Tj^, 14 1. This remarkable 
effect, in which wavefield components with zero-crossings 
at atom positions avoid photoelectron production, re- 
duces absorption in the beamsplitter by many orders of 
magnitude. Crystals XI and X2 operate at 2-beam dy- 
namical conditions, for which reflectivities of greater than 
90% are possible. For the arrangement in Fig. [T] the 
important experimental challenge is to generate three in- 
cident beams that converge focused onto a micron-sized 
volume of space at the same instant in time. "Diffraction 
before destruction" will require pulse durations of 10 fs 
or less, corresponding to a spatial pulse length of at most 
3 /zm, a technically challenging but feasible length scale 
for experimental realization. (Calculations show that 8 
keV x-ray pulses reflected from Si(lll) are stretched by 
about 4 fspH.) 

A monolithic integration of this arrangement may be 
possible. The arrangement in Fig. [2l with sample 
mounted on the beamsplitter, provides isochronal, but 
unfocused, optical paths to the sample. This arrange- 
ment is better suited to long exposures of continuous ra- 
diation for stationary samples. We do not provide de- 
tailed phase-space matching calculations here, but note 
that the beam divergence of one planned FEL, the Linac 
Coherent Light Source at Stanford, USA, is 1.1 x 10~ 6 
radians, which is less than a typical perfect-crystal rock- 
ing curve width of 3.4 x 10~ 5 radians. An energy spread 
of AE/E = 1.4 x 10~ 4 can be expected after monochro- 
mation at 8 kV, with a beam width of 20 /im. Estimates 
suggest that even with less than one scattered photon per 
pixel, phasing and reconstruction is possible [53, 



We conclude that a determination of the relative orien- 
tation between successive particles of unknown structure 
(each initially in a random unknown orientation with re- 
spect to the laboratory frame) may be achieved without 
the need to solve the phase problem. 

For a stream of identical molecules in random orien- 
tations, this would allow data from different molecules 
to be merged in the correct relative orientation. After 
phasing the resulting three-dimensional reciprocal-space 
data, a tomographic image can then be reconstructed. 
We have also shown that the orientation of successive ob- 
jects can be determined from autocorrelation functions, 



so that a solution of the phase problem is not required. 
The entire procedure cannot distinguish enantiomorphs. 
Stereoscopic projections might be obtained from just two 
projections. We have suggested experimental implemen- 
tations for this method for femtosecond X-ray diffrac- 
tion. This analysis applies to any penetrating particles 
(e.g., neutrons or high energy electrons), insofar as the 
scattering can be characterized by a scalar potential and 
the orientation and structure of the sample is unknown 
(unlike goinometer-based systems where both coordinate 
systems are known). This might include, for example, the 
tracking of the orientation of a single body from which 
nondestructive diffraction patterns can be obtained as 
function of time. "Proof of principle" measurements at 
optical wavelengths are currently under way. 

This work was supported by NSF IDBR 0555845 and 
CBST at UC Davis. 
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